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PREFACE 


This  Is  one  of  a  continuing  series  of  studies  aimed  at 
more  fully  exploiting  the  capabilities  of  modern  digital 
computers  for  the  solution  of  scientific  and  technological 
problems  of  Interest  to  the  Air  Force . 

In  this  Memorandvim  a  new  approach  to  the  problem  of  the 
determination  of  orbits  of  bodies  Is  described  and  some 
numerical  evidence  showing  Its  efficiency  is  adduced. 
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SUMMARY 

• 

Many  problems  In  mechanics  and  applied  mathematics  in 
general  require  the  solution  of  multipoint  nonlinear  boundary 
value  problems.  A  computational  approach  to  their  solution, 
based  on  the  idea  of  quasilinearization,  is  described  and  the 
results  of  an  application  to  the  computational  determination 
of  orbits  are  presented. 
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ORBIT  DETERMINATION  AS  A  MULTI-POINT  BOUNDARY  VALUE  PROBLEM 

AND  QUASILINEARIZATION 

1 .  INTRODUCTION 

The  basic  ideas  underlying  the  solution  of  boundary  value 
problems  via  quasilinearization  are  discussed  in  [l]  and  [2]. 
Results  of  some  numerical  experiments  involving  two— point 
boundary  value  problems  for  nonlinear  ordinary  differential 
equations  (Euler  equations)  are  available  in  [2]  and  [5], 
while  results  for  nonlinear  partial  differential  equations  are 
contained  in  [4].  The  purpose  of  this  note  is  to  describe  an 
application  to  the  computational  solution  of  some  general 
N— point  boundary  value  problems  and  to  present  an  application 
to  orbit  determ.lnation. 

2.  AN  N-POINT  BOUNDARY  VALUE  PROBLEM 

Consider  an  N— dimensional  vector  x(t)  =  (x^( t) ,Xp( t ) ,  .  .  . , 
which  is  a  solution  of  the  vector  system  of  equations 

(2.1)  X  =  f(x,t)  ,  0  <  t  <  b, 

and  in  addition  satisfies  the  N  conditions 
N 

(2.2)  2  aj(t^)xj(t^)  =  b^,  i  =  1,2, ...,Nj 

where 

(1.3)  0  <  t^  <  tp  <  •••  <  <  b. 

Let  us  assume  that  the  equations  (2.1)  and  (2.2)  possess  a 
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unlque  solution  on  the  Interval  [0,b].  Our  objective  Is  to 
provide  an  efficient  computational  algorithm  for  determining 
the  function  x(t)  on  the  Interval  [o,b]. 


SOLUTION  VIA  QUASILINEARIZATION 
(  \  ® 

Let  X  ^(t)  be  an  Initial  approximation  on  the  given 
Interval.  The  (k  +  1)— st  approximation,  x^^'*’^^(t).  Is  to  be 
obtained  from  the  k— th  approximation,  through  use  of 

the  relations 


1  =  1 , 2 , . . . , N  and  k  =  0,1,2,..., 


wy^h  the  boundary  conditions 

(3.2)  Z  ^ 


•  f 


N, 


for  k  —  0,1,2,....  The  linear  sys tern  for  x ^ k+1 ) ^  equat Ion 
(^.1),  subject  to  the  linear  conditions  of  equation  (5.2),  Is 
to  be  solved  In  the  following  manner.  We  determine  a 
particular  solution  of  equation  (5*1)  and  N  Independent 
homogeneous  solutions  via  numerical  Integration.  This  Is 
readily  accomplished  since  these  are  all  Initial  value  problems. 
Then  x^^'*’^^(t)  Is  expressed  as 


.(1^+1) 


(t)  = 


N 

Z 

J=1 


Vi-i) 


using  an  obvious  notation.  The  N  unknown  constants,  Cj, 
are  then  determined  by  solving  the  linear  algebraic  equations 
obtained  by  substituting  the  expression  In  (3.3)  In  equation 
(2.2). 

In  this  way,  the  solution  Is  effected  by  reducing  the 
original  problem  to  a  sequence  of  linear  Initial  value 
problems,  a  process  well  suited  for  machine  calculation. 

Some  aspects  of  convergence  (which  Is  quadratic  In  nature) 
for  the  case  of  two— point  boundary  value  problems  are  discussed 
In  [2].  We  hope  to  discuss  convergence  questions  In  a  future 
note . 

4.  ORBIT  DETERMINATION 

Consider  the  problem  of  determining  the  unperturbed  orbit 
of  a  heavenly  body  H  In  the  (x,y)  plane,  the  plane  of  the 
ecliptic,  on  the  basis  of  four  angular  observations  from  the 
earth . 


Fig.  1 — The  physical  situation 


Assume  that  the  units  are  normalized  so  that  the  equations  of 
motion  of  H  are 

(4.1)  X  =.  u,  u  -  —  x/(x^  + 
y  =  V,  0  =  -  y(x^  + 

At  times  i  =  1,2, 3, 4,  we  are  given  that 

(4.2)  e(t^)  =  e^, 
which  yields  the  conditions 

(^•3)  y(t^)  =  (x(t^)  -  l)tan  S^,  i  =  1,2, 3, 4. 

We  wish  to  determine  a  complete  set  of  elements  of  the  orbit 
in  the  form  of  values  for  x,  u,  y,  and  v  for  some 
convenient  time  t.  For  simplicity,  we  are  assuming  that  the 
earth  is  stationary  and  that  the. sole  force  on  the  heavenly 
body  is  the  gravitational  attraction  of  the  sun. 

5.  NUMERICAL  RESULTS 

Consider  the  case  for  which  the  observational  data  are 
( in  radians ) 

(5.1)  6(0.5)  =  0.251297,  6(1.0)  =  0.510240, 

6(1.5)  =  0.785690,  6(2.0)  =  1.076540. 

These  data  were  generated  by  assuming  that  at  time  zero  we  have 
x(0)  =  2.0,  y(0)  =  0,  u(o)  =  O,  v(o)  =  0.5.  The  situation, 
with  a  stationary  earth,  is  as  shown  in  Pig.  2.  As  an 


Sun 


Earth 


Fig-  2 — Approximate  directions  to  the 
heavenly  body  at  four  times  of 
observation  from  a  fixed  earth 


initial  approximation  to  the  trajectory,  we  assume  that  the 
observed  body  is  at  the  point  (1.0,  O)  at  time  zero  and  that 
it  moves  in  a  circular  orbit  with  a  period  equal  tj  3.664,  a 
very  poor  initial  approximation.  Then  using  a  Runge— Kutta  ' 
integration  method  and  the  procedure  sketched  above,  we  found 
the  following: 


Table  1 


THE  APPROXIMATIONS  TO  THE  DISPLACEMENTS 
AND  VELOCITIES  AT  TIME  2.5 


Initial 
Approx . 

First 
Approx . 

Second 
Approx . 

Third 

Approx. 

Fourth 
Approx . 

Precise 
Value  to 

6  Figures 

x(2.5) 

-0.455175 

1.554560 

1.169660 

1.190170 

1.195560 

1.195610 

u(2.5) 

-0.875520 

-1.658850 

-0.495101 

-0.664207 

-0.664265 

y(2.5) 

0.975790 

2.982690 

0.865184 

1.045950 

1.060540 

-v(g>5)j 

-0.4o8690 

1.807620, 

0.201206 

0.2576951 

0.247jt55u 

0.247499 

The  calculation  took  about  one  and  one— half  minutes,  no 
attempt  having  been  made  to  streamline  the  calculations.  It 
produced,  of  course,  not  merely  the  results  of  Table  1  but  an 
ephemerls  for  0  <  t  <  T  at  Intervals  of  0.01. 

Perhaps  even  more  suiTjrlslng  Is  the  fact  that  If,  as  an 
Initial  approximation  to  the  orbit,  we  consider  the  heavenly 
body  to  coincide  with  the  position  of  the  earth  over  all  time, 
we  still  obtain  the  same  rapid  convergence  to  the  correct  orbit. 

6 .  DISCUSSION 

The  method  Is  readily  extended  to  cases  where  more  than 
four  observations  are  available  through  use  of  the  method  of 
least  squares.  Furthermore,  perturbations  due  to  the  motion 
of  the  earth  and  the  presence  of  other  planets  are  readily 
Incorporated,  since  they  merely  result  In  modifications  of 
equations  (4.1).  Of  course,  closely  spaced  observations 
present  difficulties  In  the  numerical  determination  of  the 
constants  Cj  of  equation  (5.5).  These  matters  aind  others 
will  be  discussed  In  more  detail  later. 
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